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AN ANALYSIS OF AN 

ELECTROMECHANICAL ANALOGUE COMPUTING SYSTEM 

1. INTRODUCTION 



The object of this investigation was to 
conduct an experimental analysis and prepare a docu- 
mented study of the functioning of an electromechanical 
analogue computing system. It was the ourpose thereby 
to establish the frequency range over which the system 
operates effectively, to discover and evaluate any addi- 
tional time delays heretofore unknown, and to reveal 
any discrepancies in the oresent method of utilizing 
the system. 

The comoutlng system described and analyzed 
herein was developed because of the need of an instru- 
ment or device to solve a oarticular differential equa- 
tion. It was necessary to construct an electromechani- 
cal computer of the analogue type to simulate the rolling 
motion of a body in soace rotating about its longitudinal 
axis. The equation of motion can then be represented 
by a linear second order differential equation with 
constant coefficients, where the coefficients are functions 
of the dynamic characteristics of the particular test 
vehicle. 

The experimental results were obtained by 
varying the frequency of the sinusoidal input signal 
and studying the phase lag of the outout signal and the 
amplitude ratio of output to input signal. These closed 
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loop responses were obtained for each set of values 
selected for the coefficients of the differential 
equation. The interrelationships of these coeffi- 
cients with such parameters as system gain, damping 
ratio, and the undamped natural frequency will become 
apparent later in this report. 

Plots of these experimental results were 
then correlated with theoretical curves obtained from 
calculations based upon the assumed functioning of 
the system comoonents. 

The Investigation was carried out in the 
following stages: 

(a) Theoretical analysis of the system 
based upon the assumed functioning of the various 
components ; 

(b) Determination of the range of linear- 
ity of each individual component of the system and 
thereby the establishment of the linear range of 
operation of the over-all system; 

(c) Determination oi the operating charac- 
teristics of the closed loop system and correlation 
with the theoretical results; 

(d) Isolation and evaluation of any 
additional time delays unaccounted for in the orig- 
inal system design. 
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This investigation was conducted by the 
author at The Johns Hopkins University Aoplied Physics 
Laboratory, 8621 Georgia Avenue, Silver Spring, Mary- 
land from July, 1951 until April, 1952. 
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THEORY 

2 . I Basic Configuration . 

The electromechanical analogue computing 
system Is composed of the following basic components: 

(a) Servo amplifier. 

(b) Two-phase Induction servo motor. 

(c) Two-phase Induction generator. 

(d) Turntable or roll platform. 

(e) Gear train and other connections. 

These components are Interconnected as Indicated In 
Figure 1. 




Fig. 1. Functional Diagram of the Computing System. 

2 . 2 The Differential Equation to be Solved . 

The differential equation to be solved can be 
expressed by 

10 t G0 » -Hljj + T (2.2-1) 

where G = Angular displacement about the longitudinal 
or roll axis. 

ijj a Angle of control surface from neutral position. 
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I = Inertia. 

G • Coefficient of damping; l.e., due to wings. 

H - Forcing function; l.e., roll torque per unit 
control surface angle. 

T ■ Roll torque due to external causes. 

Neglecting T and considering a perfect 

"proportional servo" control system; 

ijj = K0 (2.2-2) 

where K Is a constant, but adjustable, ratio of control 
surface displacement to roll angle. 

The equation of motion now becomes 

e + Ge + Ike-o (2.2-3) 

I I 

or 

0 ♦ a6 + b9 = 0 (2.2-4) 

where the coefficients a and b become measures of the 
damping and stiffness resDectively . 

2.3 Determination of the System Transfer Function by 
Dynamic Analysis . 

A very brief functional description of the 
diagram in Figure 1 is as follows: The input signal (ijl) 

is generated by the movement of a control surface. It 
Is then amplified (u$) and drives the induction motor. 
The motor torque (ukljj) through an 18:1 gear ratio (n), 
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drlvea the turntable. The generator, mounted on the 
output shaft of the motor, is a damping device with an 
out out proportional to the angular velocity of the turn- 
table. Its output (f©) provides negative feedback to 
the amplifier. 

If there were no damping, the driving torque 
(uknlf) would equal the inertia retarding torque (I©) 
and the equation of motion would be 

uknty = 10 (2.3-1) 

where u is the amplifier gain, k is the motor gain, n is 
the gear ratio, ijf is the Input signal, I Is the inertia 
of the turntable and 0 is the angular displacement of 
the table. 

If now the damping effect of the generator is 
added, the equation of motion becomes 

10 = ukn (ijl - f 9 ) (2.3-2) 

Therefore , 

10 ♦ uknf© = uknljr (2.3-3) 

Dividing by I and simplifying gives 

0 ♦ a© = b$ (2.3-4) 

where a = uknf/l and b = ukn/l. It is now apparent 
that if Equation 2.3-4 is to correspond to Equation 2.2-4 
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then the coefficients a and b are functions of the 
dynamic characteristics of the particular vehicle under 
test . 

Employing the Laplace Transform (Ref. 4)^ and 
indicating the Laplace operator by "s" gives 



(s^t as) «-(s) * bljl (s ) (2.3-5) 

Thu3 , 

Q(a ) _ b b _ b/a 

f(s) s 1 * as ' s(s + a) " s (l/a s + 1) 

K 

s (Ts * 1) (2.3-6) 



This is the open loop transfer function of the system 
when the computer is one comoonent of a closed loop as 
in Figure 2. 




Fig. 2. Computer as one Component of a Closed Loop. 



References are listed in the bibliography of this 

report # 






- 8 - 



2.4 Determination of the Overall System Transfer 
Function by Servo Analysis . 

The overall transfer function of the system 
can also be found by using block diagrams and the 
Individual transfer functions. 




Fig. 3* Block Diagram of the Computer Components. 

Figure 3 is a block diagram of the components. 
If there is assumed proper design of the gear ratio 
through rigid shafting, and thus no adaltional time 
lag in driving the turntable, the diagram can be sim- 
plified to that of Figure 4 without changing the form 
of the system transfer function. 




Fig. 4. Simplified Block Diagram. 
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It Is now necessary to determine the individual 
transfer functions of the amplifier, the motor, and the 
generator. 

Amplifier . The amplifier is assumed to be 
linear, within specified limits, with negligible time 
delay. Thus, its transfer function, G-j^s), is simply K a . 

Kotor . There has been much controversy over 
the form of the transfer function of a two-phase 
Induction servo motor. The usual analysis (Ref. 1), 
based on steady state speed torque characteristics, 
results in a Laplace transfer function of the form 



G2 (s ) 



*m 

9 (TS + 1) 



(2.4-1) 



where, admittedly, -the effect of the Inductance and 
resistance of the windings has been neglected and the 
time constant is purely mechanical or a function of 
Inertia and damping. Recently there have been indica- 
tions (Ref. 5) that the transfer function might be more 
correctly expressed as 



Go (s ) = 



K. 



m 



s (T^s + 1) (T p s + 1) 



(2.4-2) 



where the second time constant may be considered an 
electrical delay equivalent to the ratio of inductance 
to resistance of the rotor. Initially it will be assumed 
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that Equation 2.4-1 Is the correct expression for the 
motor transfer function. 

Generator . Since the output of the generator 
Is proportional to the angular velocity of the turn- 
table , 

© g = K g 0 (2.4-3) 

or 

© g (s) = K g 8e(s) (2.4-4) 

Therefore, the feedback transfer function Is 

H(s) = ^S(s) . K g s (2.4-5) 



Figure 4 can now be redrawn as shown in 

Figure 5. 




Fig. 5. Block Diagram with the Associated 
Transfer Functions of the Various Components. 



From Figure 4 the overall system transfer 
function can be expressed as (Refs. 1,2, and 3) 

G 1 (s)G 2 (s) 



G(s) = 



1 + G 1 (s)G 2 (s)H(s) 



(2.4-6) 



Substituting for the Individual transfer functions gives 
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G(s) 



or, finally, 



G(s) = 



K a*m 

s (Ts + 1) 

1 + K a K mV 

a (Ts + 1) 

“a^m 

s (Ta ■» 1) + K a K m K g s 
^a^ra 

8 [Ts . (1 t ^Kgj] 



K 1 



s (T ' a ♦ 1) 



(2.4-7) 



(2.4-8) 



(2.4-9) 



(2.4-10) 



This form of the overall system transfer 
function, determined by servo analysis, Is the same 
as that determined by the dynamic analysis in Section 2.3. 

If, however, the transfer function of the 
motor were actually 



G 0 (s) = 



K 



m 



s (T x s + 1) (T 2 s + 1) 



(2.4-2) 



then from Equation 2.4-3 



K a K m 



G (s ) = 



s(T 1 s ♦ 1) (T 2 s ♦ 1) 



1 + w 



s(T lS * 1 )(T 2 s + 1) 



(2.4-11) 



KaKm 



s(T l9 + 1 )(T 2 s ♦ l) + K a K^K g s 



(2.4-12) 
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or 



G(s) = 






[Tx^bS (T x + T 2 )s + (1 + K a K m K g jj 



(2.4-13) 



K a K m 



s (s + c ) (s + d) 



(2.4-14) 



Where -c and -d are the roots of the quadratic. 



Finally, 



G(s) = 



K' 



s(T{s «t 1) (T^s + 1) 



(2.4-15) 



2.5 Theoretical Development of the Expressions for 
the Amplitude Ratio and the Phase Shift for a Second 
Order System . 



If G(s) q 



9(a) 

ijj(s) 



K 

s(Ts t 1) 



( 2 . 3 - 6 ) 



is the open loop transfer function of the system, then 
the closed loop transfer function can he expressed by 
(Ref. 1, page 141) 



Thus, 



G(s) 



G(s) 



0(°>o 

1 + G (s ) 0 

fo. _ K/s(Ts + 1) 
©i 8 1 + K/s(Ts t 1) 

K 

Ts + s + K 
1 

T 1 
K S + K S + 1 



(2.5-1) 



(2.5-2) 

(2.5-3) 



(2.5-4) 
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If and © Q , the Input and output signals, 
are sinusoidally varying quantities, the Laplace opera- 
tor "s" may be replaced by Jw (Ref. 1, page 96), where 
J is the complex operator and to is the angular sinu- 
soidal input variation. 

Therefore, 

G(Ju)) = 1 (2.5-5) 

+ £(J«) + 1 



1 

(1 -<o‘T/K) + J(^/K) 

But 

G(Jw) = ~(Jw) 
e i 



( 2 . 5 - 6 ) 



Therefore , 



;—(**» 

e i 



i 

/(I -u^t/K) 4 - (tj/K) 1 



(2.5-7) 



and if the phase between the output and input signal 
is 0 , then i- -i 



0 



arc tan 



< f/K 

1 -ufT/K 



( 2 . 5 - 8 ) 



= - arc tan 




(2.5-9) 



Where K = b/a, T = l/a and a and b are the constant 
coefficients of the equation of motion. 

Equations 2.5-7 and 2.5-9 provide a means of 
theoretically predicting the amplitude ratio of output 
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to Input signal and the phase angle by which the 
output lags the input. 

2,6 The Interrelationships of the Constant Coefficients 
and the System Parameters of Damping Ratio and the 
Undamped Natural Frequency . 

If the damping ratio, , and the undamped 
natural frequency, cu n , are to be expressed in terms of 
the constant coefficients a and b, the interrelation- 
ships existing among them must be determined. 

From Section 2.2 the equation of motion is 
a second-order equation of the form 

ID 1 * GD * Hk = 0 (2.6-1) 

or 

D 1 * aD + b = 0 (2.6-2) 

where D is the operator for differentiation with 
respect to time. This equation has two roots R^ and R 2 : 

R x and Rg = -G t |/G - 4lHk (2.6-3) 

21 

-a ^ /a - 4b 
2 

For all stable systems, that is, those whose oscilla- 
tions always attenuate to zero, the roots R^ and Rg 
may be (1) conjugate complex with negative real parts. 
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(2) equal negative real, or (3) unequal negative real, 
depending upon the numerical values of I, G, and Hk. 

Two parameters, namely, the undamped natural 
frequency, 6 J n , and the damping ratio, f f may be chosen 
so that the roots and Rg may be expressed In terms 
of two quantities Instead of three, as in Equation 2.6-3. 
One of these quantities, has the dimensions of time; 

the other one, f , Is dimensionless. Thus, if desired, 
a study of system behavior in dimensionless form can 
then be made. 

The quantity is given by the characteristic 
Equation 2.6-1 when the damping is zero, that is, 

IDS Hk » 0 (2.6-4) 

giving roots 

R x and R g - * J|/ Hk/l * k J^b" 

= ± J60 n (2.6-3) 

therefore, U) n = . The parameter f or damping ratio 

is the ratio of actual, damping to critical damping. 
Critical damping occurs when the radical of Equation 
2.6-1 is zero or when the two roots R-j_ and Rg are equal. 
The critical value of damping, f c , which can cause this 
condition is 



f c = 2J/±TET 



( 2 . 6 - 6 ) 
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The damping ratio is then 

j 3 _ actual damping m G 

f c 2|/HkI 

- a 

2|/b~ (2.6-7) 

2,1 Calculated Results from the Expressions for the 
Amplitude Ratio and the Phase Shift Versus the Input 
Frequency . 



It was desirable to be able to obtain the 
amplitude ratio, A, and the phase lag, 0 , as functions 
of the input frequency in cycles per second, for each 
set of values given to the constant coefficients a 
and b. 

A suitable range of representative values of 
a and b was then chosen. Values of a and b correspond- 
ing to damping ratios of 0.17, 0.3, 0.3 and 1.0 were 
selected and the phase and amplitude curves were calcu- 
lated for angular frequencies ranging from 0.5 through 
10.0 cycles per second. 

The 'calculated phase-lag data are tabulated 
in Table I while the amplitude ratio data are recorded 
in Table II. These results are plotted in Figures 6 
and 7» showing, respectively, the phase lag and ampli- 
tude ratio versus input angular frequency. 



-17- 



Table I . Calculated Phase-Lag Data for a 
Second Order System , 



Frequency Phase lag, 0, In degrees 



cycles/second 


i s 0.17 


±1 0-30 , 


f- 0.50 i 


= 1.00 


0.2 


- 


— _ 


..... 


11.4 


0.5 


5.2 


9.1 


14.9 


28.1 


1.0 


12.7 


21.8 


33.7 


53.0 


1.5 


30.1 


45.8 


59.6 




2.0 


90.0 


90.0 


90.0 


90.0 


2.5 


143.0 


127.0 


114.3 





3.0 


157.8 


144.2 


130.1 


112.7 


4.0 


167.2 


158.1 


146.2 


126.8 


5.0 


170.9 


164.1 


154.4 


136.5 


7.0 


174.0 


169.5 


162.8 


148.2 


10.0 


176.0 


172.9 


168.3 


157.4 


Table II. 


Calculated 


Amplitude 


Ratio Data 


for 



a Second Order System . 



Frequency 

in 




Amplitude 


ratio, A 




cycles/ second 


0-17 


/-0-30 . 


ts 

0 

• 

01 

0 

1 ? 


1.00 


0.5 


1.060 


1.050 


1.040 


0.944 


1.0 


1.300 


1.240 


1.110 


0.803 


1.5 


1.970 


1.550 


1.150 


0.641 


2.0 




1.650 


1.100 


0.502 


2.5 


1.420 


1.060 


0.730 




3.0 


0.740 


O.65O 


0.510 


0.309 


4.0 


0.340 


0.310 


0.280 


0.201 


5.0 


0.190 


0.180 


0.170 


0.138 


6.0 




0.120 




0.100 


7.0 


0.090 


0.087 


0.084 


0.075 


9.0 








0.047 


10.0 


0.042 


0,041 


0.040 


0.038 


1.97* 


2.980 








1 Q1* 




1 fi 70 










Ufw 






1.85* 






1.170 





* Denotes the frequency at which the corresponding 
amplitude ratio is a maximum. 

= 1 
p zWrT- 



(Ref. 1, page 107) 



PHASE ANGLE , <p (d e3^ees j 



-> nii-Ui.^1 Ithniu . i )v|*n •* i> to th*- iv u. ft, 
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AMPUTUDC ratio 




1 
















































r 


r-+ 
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H 


b> 
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Li 
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2.8 Development of Expressions for the Amplitude 
Ratio and Phase Shift for a Third Order System . 



If the open loop transfer function for a 
third order system Is given in the form of 2.4-15» 



G(s) 0 



K 

sO^s + 1 )(T 2 s + 1) 



( 2 . 8 - 1 ) 



then the closed loop transfer function Is given by 

G(.) = 3( ?1 9 ^ H L (2.8-2) 

1 + K 



sdjs + 1 )(T 2 s + 1) 



K 



s (T^s + 1) (T 2 s + 1) + K 



(2.8-3) 



T-iTp , Ti ♦ T;> , 1 . 

i_£ 3 3 + + -s + 1 

K K K 



(2.8-4) 



If the input signal Is sinusoidal, "s" can 
be replaced by Jcj and the equation becomes 



G(jo) = 



T-.T- 3 T- + 

4-2(Ju) ♦ -i- - 
K K 






Juj_ 

K 



(2.8-5) 



+ 1 



p.-w a /K(T 1+ T 2 )J + J |yK-0J i /K(T 1 T 2 )J 



( 2 . 8 - 6 ) 



e. 



G(Jw) = -2(Ju>) 






and 
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Therefore, 



e l 



1 

V{ -^(Tx+*a)]«^- r(V 2 )] 2 



( 2 . 8 - 7 ) 



The phase angle between output and Input signals is 
then 



0 



arctan 



co - o/TiT? 

K - coMT-L+Tg) 



( 2 . 8 - 8 ) 
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3. EXPERIMENTAL PROCEDURE AND RESULTS. 

3 ♦ I Description of Components . 

3.1-a The Power Amplifier. 

A schematic of the amplifier and power supply 
Is shown In Figure 8. • The components are built into the 
body of the computer with the exception of the mixing 
circuit which is contained in the externally mounted 
case shown in Figure 11. 

The mixing circuit is a conventional two- 
channel mixer which adds the output of the turntable, 
from a #1 Control Transformer, to the negative feedback 
of the generator; the algebraic sum feeding into the 
inverter network. 

The phase inverter receives a single-ended 
input and derives voltages for the push-pull grids that 
are 180 degrees out of phase over the complete frequency 
range of the amplifier. 

The output tubes are in push-pull to obtain 
increased power with less distortion; the even harmonics 
cancel. The power out is proportional to the number of 
tubes used. Thus, eight 6 l6's are used to give suffi- 
cient output power. 

This unit was designed and constructed at 
the Applied Physics Laboratory. 
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3.1-b The Motor. 

An Arma #6 Induction Motor (Reference 6) is 
a two-phase two-pole motor having a thin shell cup 
rotor and a stationary core, and is designed to oper- 
ate on 60 cycle alternating current. There are two 
field windings which are in space quadrature: one 

the main field, and the other the control field. A 
constant 60 cycle voltage of 115 volts is applied to 
the main field and a variable 60 cycle voltage, which 
is in phase quadrature to the main field voltage, is 
applied to the control field. Reversal of rotation is 
is accomplished by reversal of direction of the volt- 
age apolied to the control field. 

The applied voltages are rarely balanced, for 
to achieve control of torque and direction of rotation, 
the control-winding voltage may be varied from rated 
value in one sense through zero to rated value in the 
opposite sense while the voltage on the other winding 
is held constant. 

When it is necessary to have a fast acting 
yet stable closed-cycle control system, it is highly 
desirable that the torque-to- inertia ratio of the servo- 
motor have a large value. For this reason the rotor 
consists only of a thin hollow cup of light conduct- 
ing material in order to reduce inertia. The magnetic 
path is completed by means of a stationary cylinder of 
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iron placed inside the rotor. There are thus two clear- 
ance air gaos. Motors made in this manner are frequently 
called "drag cup" motors. 

This type of construction means that the 
exciting current as compared with that of motors of 
conventional design, is relatively large. Closed-cycle 
control systems normally require a servomotor with 
characteristics typical of high- rotor resistance. High- 
rotor resistance may combine with high exciting current 
to produce comditions which prevent the motor from devel- 
oping any running torque when operating single phase 
regardless of speed. This is a very desirable charac- 
teristic for many control applications although it is 
an intolerable condition in split-phase or capacitor 
motors . 

The dimensions of the Arma #6 Motor are given 
in Figure 9 while its relative size and location within 
the computer system is photographed in Figure 11. 

3.1-c The Feedback Generator. 

In servomechanisms and electromechanical 
computing systems it is often necessary to generate an 
A. C. voltage which is proportional to the speed of an 
input shaft. The ideal machine for such applications 
1 s one which generates a voltage proportional to the 
rotor speed without the introduction of phase shift, 
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The -A rrnQ 6 Induction Motor 
Fig. 9 
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harmonics or quadrature voltage. It should also be 
small in size, require little excitation power, have 
low inertia and a minimum of friction. 

Such a device is the Arma ft 3 Induction 
Generator having a stationary core and a thin shell 
cup rotor. There are two fields, main and control, 
whose axes are in space quadrature. A voltage is in- 
duced in the control field proportional to the speed 
of the cup. 

The generator dimensions are given in Figure 
10 while relative size and location are photographed 
in Figure 11. 



3.1-d The Assembled Computer. 

The entire system is housed in a device about 
40 inches wide and 36 inches high. It is portable and 
completely self contained (less recorders). Figure 11 
is a photograph of the computer. Since no pick-off 
devices were mounted during the analysis, the turntable 
has been replaced by a rotating shaft which drives the 
arm of the recording potentiometer. The control panel 
contains dials for setting the "damping" and "stiffness" 
of the system. The use of these and the various other 
dials and switches will become evident as the calibra- 
tion procedure is described. 
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The Arma 5 Induction Generator 

Fig. 10 
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3.2 Computer Calibration and Operation . 

Perhaps better understanding of the operating 
principles of the computer can be obtained by means of 
a typical illustration. 

Figure 12 is an amplified version of Figure 2 
and shows a closed loop system composed of the computer 
and a test servo. The test servo Is only an example 
and could be of many configurations to perform one of 
several functions. Possible uses are, for example; the 
simulation of a control system to stabilize a gun plat- 
form aboard ship or plane, roll stabilization of a guided 
bomb or rocket, and even the roll stabilization of a 
submarine, although In actuality other methods are 
employed. 

Perhaps one of the paramount advantages of 
this system should be emphasized here. The presence 
of the turntable or roll platform, upon which the 
actual sensing devices of a test servo-loop may be 
mounted. Is of extreme Importance. This allows exact 
duplication and minimizes the amount of analytical 
simulation necessary. If It Is desirable to study the 
effects of external disturbance torques on a particular 
test vehicle, the gyroscopes from the vehicle can be 
mounted on the turntable and the disturbance torque 
introduced into the system by merely tilting the table 
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Figure 12 . Computer and Test Servo Loot ), 

and suddenly releasing it. In this manner the computer 
can be used to study the response, in roll, to any 
external disturbance by introducing a deflection as a 
step function. 

Returning to Figure 12, the turntable carries 
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a free gyro and a rate gyro. An autosyn on the free 
gyro delivers a signal voltage proportional to the angle 
of rotation (0) of the table; and the rate gyro, through 
a pick-off, supplies a voltage in series with the former 
but of opposite phase, for damping purposes. 

Upon deflecting the platform from its neutral 
position, a signal proportional to the roll angle (0) is 
fed into the test servo; which, by means of its hydrau- 
lic piston, turns the flap through the angle ty. A 
potentiometer operated by this flap establishes a 
signal voltage proportional to ty. This voltage is fed 
to the computer amplifier, the output of which powers 
the servo motor which in turn drives the turntable and 
causes it to return to its original position. The 
system thus is a computer designed to solve the dynamic 
equation. Damping of the computer servo is effected 
by the Arma generator on the shaft of the motor, its 
output being in series with the amplifier input. 

To set up the computer for operation, it is 
first necessary to introduce values of a and b as 
determined from dynamic studies of the test vehicle 
being used. This process calibrates the system for 
operation. 

With the "calibrate-operate" switch in the 
"calibrate" position (disconnecting the test servo), 
the damping coefficient, a, is made zero (or practically 
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so) by turning the computer damping control to zero 
(adjusting Variac D). Equation 2.2-4 then becomes 

© t b0 = 0 (3.2-1) 

from which the undamped natural frequency was found to 
be Vb~ (Section 2.6), or 

= l/b/2TT cycles per second. (3.2-2) 

If then, with the damping control set at zero, the 
turntable is displaced and allowed to oscillate, its 
slightly damped frequency will be practically equal to 
the undamped frequency given by Equation 3.2-2. To 
calibrate, therefore, it is merely necessary to adjust 
the stiffness control (Variac C) until the frequency of 
oscillation, as measured on a Brush recorder, is equal 
to the value given by Equation 3.2-2. More specifically, 
a potentiometer, whose arm is driven by the turntable 
shaft, is connected to a Brush magnetic oscillograph 
which orovides an amplitude versus time trace from which 
the "natural" frequency of the computer loop may be 
determined. The voltage across potentiometer A, driven 
by the turntable, is adjusted by means of Variac C. 

This adjustment varies the voltage fed into the amplifier 
for a given angular displacement of the table; and 
hence, in effect, changes the stiffness factor, b, 
of the system. 
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To calibrate the damping, or set in the 
correct value of a, a 9teady-state roll of the turntable 
must be oroduced. To do so, the comouter servo loop 
must be broken and a signal introduced corresponding 
to a constant flaD displacement, If. This constant signal 
causes the drive motor to run at a constant speed. 

It is apparent from Equation 2.3-4- 

9 ♦ a© - bty 

that if is given some fixed displacement, l(l 0 (say 5° ) » 
the system rapidly acquires a steady roll velocity, 
whose magnitude is given by 

0 s blf 0 /a radians per second. (3.2-3) 

Since b is known and has been set in, and since is 
a known constant quantity, it is a simple matter to 
compute the value of velocity which will occur when 
the desired value of a i9 substituted in Equation 3.2-3. 
The damping calibration procedure, therefore, is to 
adjust the damping of the system, by means of the damping 
control (Variac D), to obtain a Brush recording indicat- 
ing the required value of 0. This procedure inserts the 
required damping factor, a, into the system. It should 
be noted that Equation 3.2-3 gives the velocity in 
radians per second and not in revolutions per second. 

In order to supply a constant angle of to the system, 
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the computer has a displacement signal switch (Fig. 11); 
which, when turned to either the left or the right 
position, applies an effective constant flap displace- 
ment of 5 degrees to the system. Thus, the value of ljJ o 
to use in Equation 3.2-3 is 5TT/180. When the above 
steps have been carried through, the computer is 
theoretically calibrated and ready for use. 

In the interest of learning more about the 
device, a method for checking the accuracy of calibra- 
tion will now be discussed. To do this, a recording of 
the damped oscillation or "transient response curve" 
following an initial displacement of the turntable may 
be compared with the corresponding theoretical "transient 
response curve" calculated from Equation 2.2-4. The 
solution of this equation is derived in the appendix 
and is given by 

0 = — ^2- e at//2 cos (cj t - f) ) (3.2-4) 

cos r 

where cj = /b - ( a/2 ) x 

P - arc tan a/2 

A z initial displacement of table in radians. 

°o - 

Having now a theoretical expression for the 
response, it is possible to experimentally measure the 
response by displacing the turntable, releasing it, and 
recording the motion by means of a Brush recording 
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oscillograph. Comparison of the two curves will give 
an Indication of the accuracy of calibration. 

In plotting the response curves, it is not 
necessary to plot whole curves, but merely the exponen- 
tial envelope. If Equation 3.2-4 is examined, it will 
be seen that whenever jcos (go t-/3)J = 1 , 0 falls on the 

curve 

©o -at/2 

1 e 

cos £> 

which is the equation of the envelope. These points 
occur in time whenever (cot - /3 ) = nn, or t = (n ti-(5)/o, 

where n = 0,1, 2, 3, etc. It will then be necessary, 

for the present purpose, to plot only those values of 
0 which occur when t has the above values. Experiment- 
ally, these points correspond to values of t occurring 
mid-way between successive intercepts on the time axis, 
or the peak values. 

Further simplification results when it is 
noticed that, assuming Jcos(cot -£)| = 1 (which is 

true at the points in question), the logarithm of © 
is a linear function of time. Thus, if the response 
curves are plotted, as mentioned above, on semi-log 
paper, with 9 along the log axis and time along the 
linear axis, the curves should be straight lines. 
Therefore, when the peak values of the recorded 
transient response curve of the computer are plotted 
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on the same graoh sheet with the straight "exponential 
decay" line, a direct comparison between the theoretical 
and the actual transient responses can be obtained. 

Figure 13 is a series of calibration traces 
for various values of a and b. The corresponding 
accuracy check-plots are shown in Figure 14 with the 
solid lines representing the calculated exponential 
decay, while the plotted points represent the magnitudes 
of the recorded peak values. 

A study of Figure 14 indicates that the 
computer was accurately calibrated for these particular 
values. It will be noted that the experimental values 
begin to "fall off" when the peak amplitudes become 
small or the transient has nearly died away* It is 
believed that this non-linear effect is due to friction 
and is not appreciable until the amplitude of oscilla- 
tion has decreased to less than one-third of the initial 
displacement. 
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Fig. 13. Calibration Traces for Various Values of 
a and b. 
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3.3 Establishing the Linear Limits . 

The following sequence of investigations was 
conducted to establish the linear limits of the various 
components and thus insure that all subsequent opera- 
tions were performed within the linear range of the 
system. 

3.3- a Amplifier. 

v 

To verify the assumption that the amplifier 
is a linear component, the output and Input voltages 
were recorded as functions of the angle of roll of 
the turntable. These data are tabulated In Table III. 
The results are clearly shown in Figure 15 and indicate 
linearity of amplifier output up to 11 or 12 degrees of 
table roll. Consequently, the roll angle of the table 
was hereafter confined to this maximum. 

3.3- b Motor. 

From typical motor performance curves for 
the Arma #6 Induction Motor, the stall torque is found 
to be directly proportional to the control field volt- 
age. A graph of this relationship Is taken from 
Reference 6 and reproduced in Figure 16. This property 
also has been verified at the Applied Physics Labora- 
tory. 
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Table III. Amplifier Input and Output vs. Turntable 



Roll Angle . 

Roll Angle 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 



Input Voltage 
0.030 
0.057 
0.084 
0.105 
0.132 
O.158 
0.184 
0.210 
0.235 
0.260 
0.288 
0.310 
0.335 
0.360 
0.385 



Output Voltage 

8.0 

17.0 

25.0 

33.0 

40.0 

49.0 

57.0 

65.0 

72.0 

80.0 

88.0 

95.0 

102.0 

106.0 




1 
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If the speed-voltage-torque characteristics 



of the motor were Ideal, they would be of the form 
Illustrated In Figure 17. The torque would then be 




Fig. 17 Ideal Speed- Voltage-Torque Relationships 

for the Motor 

prooortlonal to the control field voltage but independent 
of speed. Therefore, it is necessary to investigate the 
speed-torque relation of this motor. This has been done 
extensively in Reference 9 for the Arma #6 motor in a 
similar circuit. The actual speed-torque relation was 
found to be that of Figure 18. Thus, the torque is 
seen to be comoaratively independent of speed at 
low velocities and becomes more and more a function of 
speed as the velocity increases. If the speed is limited 
to one-half synchronous speed (1600 rpm) in either 
direction, the torque produced is nearly a constant. 

The following analysis will reveal the limiting effect 
of restraining the speed. 
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Converting the speed in revolutions per 
minute to degrees per second .results in 

l600rom x 360deg/rev - 533 deg/sec (3.3-1) 
60 sec/mln x 18 



where 18:1 is the gear ratio between motor and 
turntable. 

The maximum angular displacement of the 
turntable occurs at the resonant frequency which is 
aporoxlmately two cycles oer second for the theoreti- 
cal curves of Figure 7* 




From the sine wave above, it is apparent 
that the maximum velocity is equal to the peak amolitude 
multiplied by the angular velocity, or 



# 

Amplitude = A = — = 333 = 333 

w 2TTf 4TT 



(3.3-2) 



= 42.5 degrees 



Obviously, other factors of linearity will 
limit the value of table rotation below this value. 
Thus, limiting the speed of the motor to one-half 
synchronous speed will still permit sufficient angular 
displacement of the turntable. 
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3.3- c Generator. 

The assumption having been made that the 
output voltage was proportional to the input shaft rate, 
it was necessary to verify the linearity of output 
voltage versus turntable velocity for the generator. 

This was done by recording the table velocity, in cycles 
per second, on a Brush Recorder while reading the corre- 
sponding output voltage of the generator. This was done 
for a range of damping ratios from 0.1 to 1.0. The 
data are tabulated in Table IV; while the velocity, 
converted to degrees per second, is plotted versus the 
output voltage in Figure 19. 

These results establish the linear relation 
between generator feedback voltage and the roll velocity 
of the table. 

3. 3- d Recording Potentiometer. 

The linearity of the recording potentiometer 
was also verified by measuring the output signal 
voltages corresponding to various angles of roll of the 
table. 

These data are recorded in Table V and plotted 
in Figure 20. 

All subsequent utilization of the system was 
conducted within the linear range of operation as 



established herein 
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Table IV . Generator Output Voltage vs. Turntable 
Angular Velocity . 



S - 1.0 



Volts 


cos 


deg/sec 


0.100 


0.105 


38.0 


0.190 


0.233 


84.0 


0.300 


0.375 


135.0 


0.395 


0.490 


176.0 


0.495 


0.620 


223.0 


0.600 


0.750 


270.0 


0.695 


0.875 


315.0 


0.795 


0.995 


358.0 


0.895 


1.143 


412.0 


0.995 


1 .250 


45O.O 


1.090 


1.375 


495.0 


1.200 


1.500 


540.0 



Volts 


J = 0.3 
cps 


deg/sec 


0.050 


0.137 


49.0 


0.100 


0.300 


108.0 


0.200 


0.625 


225.0 


0.295 


0.935 


337.0 


0.395 


1.250 


450.0 


0.495 


1.600 


576.0 


0.600 


1.930 


694.0 



i « 0.5 



Volts 


COS 


deg/ sec 


0.218 


0.110 


78.0 


0.357 


0.195 


130.0 


0.532 


0.285 


192.0 


0.667 


0.355 


240.0 


0.800 


0.425 


288.0 


1.020 


0.545 


368.0 


1.286 


0.675 


464.0 


1.470 


0.770 


529.0 


1.670 


0.870 


600.0 



Volts 


:f - O.l' 

CPS 


deg/ sec 


0.050 


0.400 


144.0 


0.075 


0.625 


225.0 


0.100 


0.855 


308.0 


0.125 


1.110 


400.0 


0.150 


1.330 


478.0 


0.180 


1.600 


575.0 



( 



- 51 - 



Table V . Output Voltage of Recording Potentiometer vs . 
Roll Angle of Turntable * 



Output voltage 
0.09 

0.15 

0.23 

0.30 

0.38 

0.46 

0.53 

0.60 

0.69 

0.74 

0.83 

0.89 

0.99 

1.04 

1.12 

1.23 

1.32 

1.37 

1.45 

1.53 

1.62 



Koll Angle 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 



21 



OUTPUT. VOLTAGE. 
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3.3-© Frictional Non-Linearity for Small Angular 
Velocities . 



From the slooe of the curve in Figure 20, the 
value of output voltage oer degree of table roll can 
be determined ( 0,075 volt per degree). 

Values of output voltages can now be converted 
to degrees of roll of the table, while the angular 
velocity in degrees per second can be calculated from 
Brush recordings. The steady-state roll velocity- is 
given by Equation 3.2-3 as 

0 = bty Q /a radians per second 

where t)l 0 is a fixed displacement of the table. 

If the angular velocity is plotted against 
table displacement, a straight line should result from 
the above linear relationship. Figure 21 shows the 
actual curve, which is linear except for small dis- 
placements at low frequencies; the velocity being a 
function of both frequency and the amplitude of dis- 
placement. From a study of Figure 21, the presence of 
frictional dead-space is apparent. This simply means 
that if experimental data are taken for small angular 
velocities of the turntable, the frictional effects 
will be appreciable and will create phase lag in 
excess to that predicted. This effect can be explained 
as follows.- 
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The slope of the theoretical curve Is 




a/b &, the damping factor. 



At low velocities, the slope of the vector from the 
origin to the actual curve Is greater, indicating the 
presence of the added frictional damping. In other 
words, due to the frictional non-linearity a greater 
displacement Is necessary to obtain the same angular 
velocity. 
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Table VI . Angular Velocity vs. Table Displacement . 



Angular Velocity 



Displacement 



Cycles 

per 

second 



Degrees 

per 

second Volts Decrees 



0.067 

0.159 

0.237 

0.345 

0.526 

0.769 

0.995 



24.1 


0.095 


1.27 


55.4 


0.165 


2.20 


85.2 


0.230 


3.07 


124.2 


0.315 


4.20 


189.0 


0.475 


6.31 


277.0 


0.685 


9.13 


358.0 


0.875 


11.61 


396.0 


0.965 


12.86 



1.100 



“ 5*6 ” 








- 
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3.4 The Closed-Loop System and Measuring Apparatus . 

The computer, connected as a closed loop, 
along with the necessary apparatus for measuring the 
system output to Input amplitude ratio and the phase 
lag. Is shown schematically in Figure 22. A photograph 
of this experimental set-up is shown in Figure 23. 

3.4- a The Computer as a Closed- Loop. 

The computer amplifier is connected to a 
110-volt 60-cycle per second power supply. The ampli- 
fier output signal feeds into the control field of the 
motor, while a constant 60-cycle per second voltage of 
110 volts is applied to the main field. The motor shaft, 
through an 18:1 gear train, drives the turntable. The 
#1 control transformer, mounted on the turntable shaft, 
returns a signal to the mixing circuit proportional to 
the rotation of the turntable. The mixer receives a 
second signal from the feedback generator which is pro- 
portional to the speed of rotation of the table. The 
mixer output is subtracted from the system input and 
the resulting signal feeds back into the computer 
amplifier. 

3.4- b The Signal Generating and Phase-Lag 
Measuring System. 



The signal generating and phase-lag measuring 
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system consists of a pair of #1 control transformers, 
two pair of autosyns, two similar d-c amplifiers, two 
similar demodulators, and a cathode-ray oscilloscope 
modified for d-c signals, ^his system, indicated in 
Figure 22, is shown in detail In Figure 24. 

Autosyns 1 and 2 are excited by 26 volts at 
8000 cycles per second from an audio signal generator. 
Control transformer 1 is excited by a 60-cycle per 
second signal from the amplifier. The rotors of control 
transformer 1 and autosyn 1 are mechanically coupled 
and can be driven by a variable-speed motor. The body 
of control transformer 1 can be rotated by means of a 
phase-shift knob, and since the case of autosyn 1 is 
fixed, this causes a phase shift between them. A 
pointer on the phase-shift knob traverses the face of 
a circular dial calibrated in degrees. 

Yfith the phase-shift knob adjusted to zero 
on the dial, the voltage across control transformer 2 
is adjusted to null by rotation of its input shaft. 

To obtain the ,, ln-phase ,, condition of autosyns 1 and 2, 
the voltmeter was connected across the autosyn pair and 
again adjusted to a voltage null. This was repeated 
for autosyns J> and 4 after setting the table to zero 
degrees of rotation. 

The modulated signal of feedback autosyn 4 is 
fed to the d-c amplifier and demodulator (2). The 
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output signals of the two d-c amplifiers and demodu- 
lators are applied to the cathode-ray oscilloscope 
modified for d-c signals. Since these two signals are 
sinusoidal, a Lissajou figure appears on the cathode- 
ray tube which is a single loop. When the two sinu- 
soidally modulated signals are in phase, the Lissajou 
figure closes to a single straight line. 

When the phase-shift knob is set to zero on 
the dial, the modulated signals from control transformer 
2 and autosyn 2 are in phase for all driving speeds. 

A phase lag occurs through the computer loop, and the 
output-modulated signal from autosyn 4 will lag the 
input signals. If the phase-shift knob of control trans- 
former 1 is rotated in such a manner as to cause the 
input signal of control transformer 2 to lead the input 
signal of autosyn 2 by an amount equal to the phase lag 
that occurs in the comouter, modulated signals from 
autosyns 2 and 4 are in phase. 

This in-phase condition is readily apparent 
when the loop of the Lissajou figure on the cathode-ray 
tube becomes a straight line. Thus, to determine the 
phase lag, it is only necessary to turn the phase-shift 
knob until the lissajou figure becomes a straight line. 
The measured phase lag is then read directly from the 



circular dial 
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3.4-c The Amplitude Measuring System. 

The amplitude measuring system consists of 
a recording ootentiometer , driven by the turntable shaft, 
a d-c supply, a zero-positioning potentiometer, a d-c 
Brush amplifier and a Brush recorder. This system is 
shown in Figure 25. 



Table shaft 



/ 

/ 

/ 




Recording 

Potentiometer 



d-c Supply 



Zero-positioning 

Potentiometer 




Brush 

d-c 

Amplifier 



Recorder 



Brush 



Fig. 25. The Amplitude Measuring System 
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The d-c supply waa connected to the recording 
potentiometer in parallel with the zero-positioning 
potentiometer. The zero-poaitionlng potentiometer was 
then adjusted ao that the potential, and thus the out out 
of the recording potentiometer, was zero when the table 
was at rest at zero degrees of rotation. The d-c Brush 
amplifier was used to amplify the output signal from 

i 

the recording potentiometer. Finally, a graph of the 
amplitude as a function of time was recorded on paper 
tape in the Brush recorder. 

The Brush recording not only supplies a plot 
of the output amplitude, but also orovides a means of 
measuring the actual input signal frequency. From the 
recorder tape, the number of cycles and the elapsed time 
are both obtained to determine the frequency in cycles 
per second. 

At very low frequencies the amplitude ratio 
of output to input is essentially unity. Therefore, the 
amplitude of the input signal was adjusted while the 
table was driven at the slowest possible speed. This 
was done by adjusting the gain of the Brush amplifier 
until the Brush recorder pen oscillated over a peak-to- 
peak distance of the desired amount. 

3.5 The Experimental Results . 

The experimentally determined amplitude-ratio 



and phase-lag data are recorded in Tables VII to X 
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These tests were conducted with the comouter calibrated 
for the same values of a and b, and thus the same damp- 
ing ratios, as were used in the theoretical calculations 
of Section 2.7. This facilitates correlation of theore- 
tical and experimental results. The phase-lag data are 
plotted versus frequency in Figures 26 and 27, and the 
amplitude ratio data are plotted versus frequency in 
Figure 28. 

The re-plotting of Figure 26 was necessary 
because of the phenomena discussed in Section 3.3-©. 

The discrepancies between the experimental and theore- 
tical phase lags at low frequencies was not believed 
to be accurate. Upon investigation it was discovered 
that if the initial displacement of the table was 
increased, the additional phase lag, due to frictional 
effects, was eliminated. The corrected data was then 
plotted in Figure 27 showing the close correlation 
expected at lower frequencies. 

Figures 27 and 28 indicate the effect of 
varying the damping ratio, at a constant resonant 
frequency. To obtain a more complete picture of the 
system performance, a sequence of tests was conducted 
at various resonant frequencies with a fixed damping 
ratio of 0.17. These data, recorded in Tables IX and 
X, are plotted in Figures 29 and 30. 
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Table VII , Experimental Phase- Lap; Data for Varying; 
damping; Ratios at Constant Resonant Frequency , 



1 = 0.17 



/ s 0.30 



Frequency 


Phase Lag 


Frequency 


Phase Lag 


in 


in 


in 


in 


cycles/second 


degrees 


cycles/second 


degrees 


0.25 


8 


0.34 


8 


0.57 


12 


0.55 


13 


0.90 


17 


0.84 


20 


1.25 


24 


1.38 


38 


1.70 


44 


1.60 


51 


1.80 


55 


1.86 


76 


1.92 


77 


2.00 


92 


2.00 


93 


2.15 


114 


2.10 


103 


2.50 


137 


2.16 


119 


3.00 


152 


2.35 


141 


3.46 


162 


2.61 


156 


3.93 


168 


3.22 


169 


5.00 


176 


3.90 


175 


O. 30 * 


5 


5.00 


180 


0 . 50 * 


9 


0 . 50 * 


5 


0 . 70 * 


14 


1 . 00 * 


13 


1 . 00 * 


21 


1 . 20 * 


18 






/ = 0 . 


50 


f = 1 . 


00 


Frequency 


Phase Lag 


Frequency 


Phase Lag 


In 


in 


in 


in 


cycles/ second 


degrees 


cycles/second 


degrees 


0.25 


8 


0.23 


13 


0.50 


16 


0.44 


25 


0.80 


26 


0.67 


36 


1.05 


36 


0.91 


49 


1.35 


51 


1.20 


62 


I .65 


68 


1.48 


73 


I .85 


84 


1.75 


82 


2.20 


102 


2.00 


93 


2.45 


115 


2.50 


106 


2.86 


132 


3.00 


118 


3.30 


144 


3.60 


128 


4.00 


155 


4.24 


138 


4.71 


162 


4.80 


146 


5.62 


169 


5.94 


155 


6.69 


176 


7.65 


166 






8.75 


172 



* Corrected data 
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Table VIII . Experimental Amplitude-Ratio Data for 
Varying Damping; Ratios at Constant Resonant Frequency . 

/ = 0.17 / = 0.30 



Frequency 




Frequency 




In 


Amplitude 


In 


Amplitude 


cycles/second 


ratio 


cycles/ second 


ratio 


0.25 


1.010 


0.34 


1.020 


0.57 


1.080 


0.55 


1.080 


0.90 


1.250 


0.82 


1 . loO 


1.25 


1.600 


1.38 


1.520 


1.70 


. 2.500 


1.60 


1.750 


1.80 


2.900 


1.86 


1.900 


1.92 


3.300 


2.00 


1.850 


2,00 


3.500 


2.15 


1.650 


2.10 


3.040 


2.50 


1.200 


2.16 


2.600 


3.00 


0.700 


2.35 


1.900 


3.46 


0.455 


2.61 


1.250 


3.93 


0.305 


3.22 


O.56O 


5.00 


0.165 


3.90 


0.335 


7.08 


0.082 


5.00 


O.170 






7.02 


0.086 






J = 0 


.50 


/ = 1 . 


00 


Frequency 




Frequency 




In 


Amplitude 


In 


Amplitude 


cycles/second 


ratio 


cycles/ second 


ratio 


0.23 


O.985 


0.25 


1.000 


0.44 


0.970 


0.50 


1.035 


0.67 


0.910 


0.80 


1.080 


0.91 


0.860 


1.05 


1.120 


1.20 


O.78O 


1.35 


1.180 


1.48 


0.705 


1.65 


1.230 


1.75 


0.620 


I.85 


1.205 


2.00 


0.540 


2.20 


1.000 


2.50 


0.420 


2.45 


0.820 


3.00 


0.320 


2.86 


0.610 


3.60 


0.240 


3.30 


0.450 


4.24 


0.180 


4.00 


0.270 


4.80 


0.145 


4.71 


0.185 


5.94 


0.095 


5.62 


0.130 


6.67 


0.080 


6.69 


0.087 


7.65 


O.O56 


8.23 


0.054 


8.75 


0.043 
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Table IX . Experimental Phase-Lap, (0) Data for a 
Constant Damping. Ratio at Various Resonant Frequencies . 



Cos 


Degrees 


Cps 


Degrees 


Cps 


Degree 


0.20 


10 


0.50 


5 


0.14 


1 


0.30 


19 


1.00 


13 


0.78 


4 


0.39 


56 


1.20 


18 


1.50 


7 


0.48 


84 


1.70 


44 


2.34 


13 


0.55 


107 


1.80 


55 


2.80 


17 


0.60 


127 


1.92 


77 


3.33 


24 


0.70 


145 


2.00 


93 


3.57 


30 


0.80 


154 


2.10 


103 


3.89 


41 


0.90 


159 


2.16 


119 


4.28 


57 


1.20 


168 


2.35 


141 


4.50 


74 






2.61 


156 


4.67 


92 






3.90 


175 


4.85 


116 






5.00 


180 


5.25 


136 










5.46 


150 










5.95 


160 










6.6 5 


171 


Table X. 


Experimental Amplitude-Ratio 


(A) Data 


for a 


Constant 


Damping 


Ratio at 


Various Resonant Frequencies. 


CpS 


A 


Cps 


A 


Cps 


A 


0.20 


1.20 


0.25 


1.010 


0.14 


1.001 


0.30 


1.55 


0.57 


1.080 


0.78 


1.010 


0.39 


2.10 


0.90 


1.250 


1.50 


1.050 


0.48 


2.80 


1.25 


1.600 


2.34 


1.300 


0.55 


2.30 


1.70 


2.500 


2.80 


1.480 


0.60 


1.60 


1.80 


2.900 


3.33 


1.800 


0.70 


0.90 


1.92 


3.300 


3.57 


2.150 


0.80 


0.55 


2.00 


3.500 


3.89 


2.660 


0.90 


0.40 


2.10 


3.040 


4.28 


3.500 


1.20 


0.18 


2.16 


2.600 


4.50 


4.200 






2.35 


1.900 


4.67 


4.330 






2.61 


1.250 


4.85 


3.700 






3.22 


O. 56 O 


5.25 


3.000 






3.90 


0.335 


5.46 


2.205 






5.00 


0.170 


5.95 


1.500 






7.02 


0.086 


6.65 


O. 8 O 5 










7.90 


0.400 



Damping ratio of 0.17. 

Resonant frequencies of 0.5, 2.0, and 4.5 respectively. 
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4 . CORRELATION OF THE EXPERIMENTAL AND THEORETICAL 
RESULTS. 

4 . 1 Experimental ver9u9 Second-Order Theory . 

The experimental phase-lag curves for the 
closed-loop system are plotted in Figures 26 and 28, 
from the data of Tables VII and IX. The phase lag in 
degrees is plotted as a function of the input signal 
frequency in cycles per second. Two sets of curves 
were obtained; one at constant resonant frequency, with 
damping ratios varying from 0.17 through 1.0, and the 
other at a constant damping ratio, with the resonant 
frequency varying from 0.5 to 4.5 cycles per second. 
Theoretical second-order system phase-lag curves, with 
these same values of j and co , are also plotted on 
Figures 26 and 28. 

At low frequencies, there is a high degree 
of correlation between the experimental and the theoret- 
ical results. However, as the frequency increases 
beyond resonance, the experimental phase shift lags, 
more and more, that predicted by the second-order system 
theory. This deviation clearly indicates the presence 
of an additional linear delay in the system, since the 
increasing phase lag appears as a function of frequency. 

The experimental and theoretical second-order 
system amplitude-ratio curves are plotted in Figures 28 
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and 30, from the data of Tables VIII and X. The 
experimental results Indicate greater over-shoot 
than expected, since the maximum values of the amplitude 
ratio are greater than oredicted. Beyond the resonant 
frequency, the negative slopes of the experimental curves 
are greater than the predicted slopes. Also, from Figure 
30 , as the resonant frequency increases, the peak values 
of the amplitude ratio Increase. All of these results 
lend credence to the existence of an additional linear 
delay, said delay being a function of frequency. 

It is difficult to quantify the effectiveness 
of the system, since the results are dependent on the 
calibration. The damping ratio, j , is a function of 
the calibration constants, a and b. Thus, as the stiff- 
ness and damping of the system are altered, the results 
vary. However, the following conclusions seem Justified 
from the results obtained: 'No appreciable error appears 
in the phase lag below 2 cycles per second. For the 
frequency range from 2 to 5 cycles per second, the error 
increases from zero to 6 or 7 per cent. Beyond this 
frequency, the error increases more rapidly. 

A . 2 Experimental versus Third-Order Theory . 

Assuming that the additional time delay of 
the system is in the form of a distinct, independent 



time constant, its magnitude can be determined readily 
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by trial and error approximations. This can be done 
by substituting values for into the third-order 
expressions developed in Section 2.8, until these 
calculated results fit the experimental curves. 
and K retain their original values, since they are 
functions of the calibration constants, a and b. 

The value of T 2 was first estimated at 3 
milliseconds and found to be too small. After two 
more trials, the value of 5 milliseconds resulted 
in a calculated point that plotted directly on the 
corresoonding experimental curve. Further calcula- 
tions led to a remarkable correlation between the 
third-order and experimental curves. This was true 
for both the phase-shift and amplitude-ratio curves. 

The calculated data are tabulated in Tables 
XI and XII and olotted with the experimental results 
in Figures 31 and 32. 



I 
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Table XI. Calculated Phase-Lap; Data for a Third 
Order gyatem . 



Frequency 

In 


Phase Lag, 


In Degrees 




Cycles/ second 


7 -0.17 


t-0.30 


f-0-50 


f *1.00 


0.2 


2.0 


3.5 


5.8 


11.4 


0.5 


5.1 


9.2 


14.9 


28.2 


1.0 


12.3 


21.5 


33.9 


53.8 


1.5 


28.2 


45.6 


61.1 


75.0 


2.5 


152.0 


135.0 


120.6 


108.0 


3.0 


167.0 


152.8 


137.7 


119.0 


4.0 


176.7 


168.0 


155.3 


135.2 


5.0 




175.0 


165.1 


146.6 


7.0 






176.6 


162.0 



Table XII . Calculated Amplitude-Ratio Data for a 
Third-Order System. 



Frequency 

In 




Amplitude 


Ratio, A 




Cycles/Second 


1 = 0.17 


i- °-3° 


i * 0.50 


II 

h-» 

• 

O 

O 


0.5 


1.065 


1.060 


1.035 


0.950 


1.0 


1.310 


1.250 


1.130 


0.820 


1.5 


2.060 


1.670 


1.210 


0.680 


2.0 


3.590 


1.860 


1.120 


0.520 


2.5 


1.490 


1.130 


0.770 


0.410 


3.0 


0.750 


0.670 


0.535 


0.320 


4.0 


0.325 


0.310 


0.285 


0.205 


5.0 


O.I 85 


0.180 


0.170 


0.140 


7.0 


0.088 


O.O 85 


0.082 


0.070 


10.0 


0.040 


0.038 


0.037 


0.035 
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5. ISOLATION OF THE ADDITIONAL DELAY. 

Discussion of Delay . 

Having established the presence of an 
additional time lag, an attempt to isolate and measure 
the delay Is believed to be incumbent upon an investi- 
gation such as this. The assumption has been made, and 
partially substantiated, that the delay is distinct and 
Independent. This indicates the already logical source 
of the additional delay--the Induction motor. 

In Section 2.4, two forms were given for the 
transfer function of an a-c , two-phase, induction motor. 
The second of these, Equation 2.4-2, contained two 
distinct delays, one mechanical and one electrical. 

The first, Equation 2.4-1, has been used in this 
analysis, until now, for two reasons: (1) The conven- 

tional texts, for the most part, have treated the 
electrical constant as negligible or have ignored it. 

(2) Utilization of the system to date has been based 
upon the system transfer function developed in Section 
2.3 » which is devoid of any electrical delay. However, 
now that a third-order system is apparent, the presence 
of an electrical delay in the motor seems very likely. 

In Reference 5* Mr. Brown substantiates this 
assumption and states the magnitude of the delay as 
equivalent to the ratio of rotor inductance to rotor 
resistance. A similar transfer function for a d-c 
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motor (Reference 1, page 128) gives the second or 
electrical delay as the ratio of the field inductance 
to resistance. Dr. Sobczyk, in his analysis in 
Reference 10, concludes that the a-c motor is effective- 
ly equivalent to a demodulating device followed by a 
d-c motor. Thus, if an analogy is drawn to the d-c 
treatment of Brown and Campbell (Reference 1, page 128), 
the following equivalent circuit results from the motor 
as connected to the amplifier. 




Fig. 33. Equivalent Circuit of the Motor and 
the Amplifier Outout Stage. 



The transfer function of such a system is 
of the form, 



6, . _ K 

E 0 " s(T m s-l) (T e s-*-l ) 



(3.1-D 



where T ffi is the mechanical delay due to the inertia 
and viscous friction, while T e is an electrical delay. 
Considering only the electrical circuit, 

/<e 3 = r p i p + Ri P + L di ? /dt 



( 5 . 1 - 2 ) 
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which is transformable to 

^uE(s) = I ( s ) [r p «• R + Ls] (5.1-3) 

Thus , 



l yU. 

E^ 3 (r P +R)+ Ls 



(5.1-4) 



Since the rotor torque ls proportional to the current, 



T / x K K 

E 18 ' (L/r P +R)s + 1 “ T e s + 1 



(5.1-5) 



Obviously, the process of calculating this 
delay would be difficult and laborious; and since the 
magnitude of the delay ls believed to be approximately 
5 milliseconds, a method of measuring it experimentally 
was devised. 

5.2 Experimental Procedure . 

To effect the measurement of the electrical 
delay it was necessary to attach a transducer element 
to the output shaft of the motor. This enabled record- 
ing of the torque-to-input transfer function of Equa- 
tion 5.1-5. 

The transducer used was a Model Gl-48-675 
of the Statham Laboratories (Reference 11). 

The Statham Transducer Element . When physical 
measurements are made it is often necessary to transduce 
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mechanical parameters to sensible electrical units. 

In force or torque studies, a displacement must actuate 
a transducing element which will yield an electrical 
signal proportional to the torque. The electrical sig- 
nal can then be transmitted to an electrical recording 
system. To satisfy these needs, Statham Laboratories 
has developed a basic displacement transducer in which 
displacements up to 0.0015 inches are caused to produce 
sensible resistance changes proportional to the dis- 
placement. 




Fig. 34. The Statham Transducer 
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The displacement transducer element Is shown 
In Figure 34. It consists of a frame, E, supporting 
a movable armature, F, by two thin cantilever plates. 

Four sets of constantan filaments, A, B, C, D, are strung 
under initial tension between the frame and armature in 
the manner shown. When the armature is displaced long- 
itudinally, two of the sets of filaments will be elong- 
ated, while the other two sets will shorten. The elong- 
ated filaments will increase in resistance, while the 
resistance of the shortened filaments will decrease. 

The change in resistance of the filaments will be pro- 
portional to their change in length. The transducer 
is so wired that the filaments are connected in a Wheat- 
stone bridge circuit, as shown in Figure 35* The resist- 
ance change of the filaments alters the electrical bal- 
ance of the bridge so that an electric current is 
caused to flow in the output circuit. 




Fig. 35. The Bridge Circuit of the Transducer 
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The cantilever plates which suoport the 
armature are mounted on the under side of the element. 

G is the linkage pin by which movement is transmitted 
to the armature. A fixed stud on the armature rides 
within a square hole in the undercarriage, so accurately 
sized that the armature can move a maximum of 0.0015 
inches in either direction. A "trimmer" resistor is 
placed in series with one of the bridge elements. It 
is adjusted after assembly of the transducer to equalize 
the resistance of the four elements of the bridge. If 
the bridge is balanced at zero output before the dis- 
placement is applied, the unbalance electrical outout 
of the bridge will bear an accurate, linear relationship 
to the displacement . 

Procedure . The transducer was attached to the 
motor shaft as shown by the photograph of Figure 36. 

The rotor torque is transmitted through the linkage 
pin and produces a displacement of the armature which 
is directly proportional to the torque. In the model 
used, a torque of 48 ounce-inches is required to displace 
the armature 0.0015 inches. 

It should be emphasized here that the maximum 
displacement of the rotor shaft is 0.0015 inches. Thus, 
the motor is essentially stationary and independent of 
the forces of inertia and viscous friction. In other 
words, if the torque equation is considered; 
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T = 19 + f0 + k© (5.2-1) 

It Is apparent that the torque from the Inertial and 
frictional terms provides a negligible contribution 
when compared to the third term. 

The amplitude and phase measuring systems 
are essentially the same as described in Section 3.4. 
The out out of the transducer element is fed through a 
d-c amplifier and demodulator to the cathode-ray 
oscilloscope, where its phase is compared with the 
phase of the input signal. The Brush recorder was 
replaced by the Miller recording system, since the 
frequencies were higher. At the higher frequencies, 
the Brush recorder contributes to the phase lag, 
while the Miller recorder has the following desired 
characteristics: The amplitude response is flat 

and the phase lag negligible over the frequency 
range used. 

The calculated and experimental data are 



tabulated in Table XIII and plotted in Figure 37 
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-Fig. 36 . T we. Mounted Transducer 
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Table XIII . Calculated and Experimental Data for a 
First-Order System . 



A = ]/\ + ( coT) 1 ^ S " arctan6jT T = °-°° 5 



CALCULATED 



f 


GJ 


A 


ft 


1 


6.28 


1.000 


1.8 


3 


18.84 


.996 


5.4 


3 


31.40 


.990 


8.9 


10 


62.80 


.955 


17.4 


15 


94.20 


.905 


25.2 


20 


125.60 


.847 


32.1 


25 


157.00 


.787 


38.1 


30 


188.40 


.730 


43.3 


AO 


251.20 


.623 


51.4 





EXPERIMENTAL 




f 


A 


ft 


2.0 


.998 


4.0 


4.5 


.995 


7.5 


8.5 


.960 


14.5 


11.0 


.940 


18.0 


14.0 


.910 


23.0 


19.0 


.850 


30.0 


23.0 


.800 


35.0 


28.0 


.750 


40.0 


30.0 


.720 


43.0 




CPS 
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5*3 Correlation of Results. 

From Figure 37 , the correlation of the 
calculated and experimental results appears to 
be very good. The physical limitations of the 
experimental apparatus allowed measurements to 
only 30 cycles per second. However, this is 
sufficient to verify the magnitude of the delay. 

For a fir3t-order system, a phase lag of 45 degrees 
occurs at the frequency, >-’ = l/T. The amplitude 
attenuates 3 decibels, or to approximately 0.7 
of its original magnitude, at this same frequency. 
When T a 0.005* ^ = 200 and, therefore, the 45 degree 
phase shift occurs at approximately 32 cycles per 



second. 
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6. CONCLUSIONS AND RECOMMENDATIONS 

As a result of this investigation, the 
author has reached the following conclusions: 

(1) The correlation between the experimental 
results and those calculated for a second-order system 
is very good for frequencies to 2 cycles oer second. 
Above this frequency, the experimental Dhase lag is 
greater than predicted, but remains less than 10 per 
cent in error out to 5 cycles per second. Above 5 
cycles per second, the discrepancy increases rapidly. 

(2) The order of the computing system is 
greater than the second, as indicated by the phase 
lag and output signal attenuation being greater than 
predicted for a second-order system. 

(3) The presence of an additional linear 
delay indicates the computer is a third-order system. 

(4) The additional delay is an "electrical" 
delay within the induction motor and has a magnitude of 
5 milliseconds. 

The following recommendations are respect- 
fully submitted: 

(1) Additional phasing capacitors should be 
added to the Arma motor circuits to insure that the 
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torque Is independent of speed over the useful speed 
range of the motor (Reference 9). 

(2) A study should be made with the objective 
of eliminating or compensating for the additional delay. 
Perhaps one approach is to increase the impedance of 
the amplifier stage and thus, decrease the ratio of 
inductance to resistance. 

(3) In order to quantify the effectiveness 
of the computer as a second-order system, additional 
sequences of frequency responses are necessary over 
varying ranges of frequency and damping ratios. Perhaos 
a worthy extension of this investigation might be the 
interpretation of these data. If it is possible to 
interpolate between such curves, then graphs or tables 
could be prepared to assign a definite index of effec- 
tiveness for a given set of operating conditions. 
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APPENDIX A: Symbols. 

Amplitude ratio of output-to-input signal 

Calibration constant measure of damping G/l 

Calibration constant measure of stiffness — H/l 

Operator for differentiation with respect to 
time 

Generator feedback 
Critical damping 
Coefficient of damping 
Open-loop transfer function 
Closed-loop transfer function 
Forcing function 

Transfer function of feedback element 

Inertia 

Plate current 

Complex operator 

Closed-loop gain constant 

Amplifier gain constant 

Motor gain constant 

Generator gain constant 

Proportionality constant 

Inductance 

Maximum amplitude ratio 
Gear ratio 
Resistance 
Roots of quadratic 
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P 

rpm 

s 

Syn 

T 

T 



L m 



u 

V c 

0 

• 

© 

© 

0 3 

9 c 

0 



Uo 

CO 



<o. 



n 



f 



r^J 



Plate resistance 
Revolutions per minute 
Laplace operator 
Synchronous 
Time constant 
Electrical time constant 
Mechanical time constant 
Amplifier gain 
Control field voltage 
Angular displacement 
Output signal velocity 
Output signal acceleration 
Input signal 
Output signal 

Phase angle between input and output signals 
Angle of control surface displacement 
Fixed control surface displacement 
Angular frequency in radians per second 
Undamped natural frequency 
Damping ratio 
Amplification factor 
Proportional to 
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appendix B: Solution of the Equation of Motion. 

The equation to be solved is 
0 + a0 + b0 s 0 

This equation must be solved subject to the condition 
of an initial displacement; l.e., at t®0, 0=0 O , and 
9=0. 

Writing the equation in symbolic form using 
the operator (D=d/dt ) gives: 

(D 2 + aD + b)9 = 0 

Therefore, 

- ~ a */a - 4b 
D 2 

When 4b > a, as is generally true, 

D = -a/2 ± i /b - (a/2 / = -a/2 ± lGJ 

where Go = /b - (a/2) z 
The solution is, therefore: 

9 = c lS ( - a/2 + lu) i Cj e ( - a/2 ' lw)t 

= cos cot + Cosin cot) 

Using the first initial condition, that 0=0 O at t=0, 
gives : Ci * 9©. 
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Similarly, the second initial condition gives; C 2 = a0 o /26J 
The complete solution then becomes: 

0 = e" a / 2 1 (0 o cos wt 4 a9 0 /2(o sin cot ) 

From the relationship, 

acos0 4 bsln0 = 1/a 1 + b 1 cos (0-tan b/a) 

the equation becomes: 

q = e " a / 2 1 /SJ7( a/^feJ^/coaCcot-^) 

\ 

where ^ = tan ^a/20J 
Thus, 

Q - e~ a/> '2 ^ e 0 /2Col/4iu 1 ta 1 cos (tot-|6) 

and since 0J= )/b- (a/2) Z 

0 = ^ 0 o l/b/ko cos(cot-P) 

If, CO = 1/b- (a/2) 2 
00 2 = b - (a/2 f 

I JL 2L 

and b = (j + (a/2) 

This relation may be represented as shown in 
the figure below: 
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From this figure, it may be seen that 
/b/u) = l/cos fb 



Thus , 

© = 0 o /cos£> e _a//2 ^ cos(oJt-p) 

where W * Vb-(a/2) x 

P - tan" 1 a/2 CO 



The damped frequency of oscillation of the 
system can be stated as follows: 



f = (0/2TT = j/b- (a/2 J^TT 
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